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Abstract 

We present a unified approach which is accurate at all temperatures and densities for 
calculating the energy loss from a stellar plasma due to the plasma process, the decay 
of photons and plasmons into neutrino pairs. To allow efficient numerical calculations, 
an analytic approximation to the dispersion equations for photons and plasmons is 
developed. It is correct to order a in the classical, degenerate, and relativistic limits 
for all momenta k and is correct at small k for all temperatures and electron densities. 
Within the same approximations, concise expressions are derived for the transverse, 
longitudinal, and axial vector components of the neutrino emissivity. 



The emission of neutrinos can be an important energy loss mechanism for very hot or 
dense stars. The collective effects of the stellar plasma can significantly alter the production 
rate of neutrinos. The most dramatic example is the "plasma process" , the decay of photons 
and plasmons into neutrino pairs, a process that owes its very existence to plasma effects. 
It was pointed out by Adams, Ruderman, and Woo |T[ in 1963 that the plasma process 
could be the dominant energy loss mechanism for very hot and dense stars. It was recently 
shown that the formulae for the plasma process that have been used in all previous work 
are inaccurate at relativistic temperatures and electron densities [[|, underestimating the 
emissivity by a factor as large as 3.185. The purpose of this paper is to provide a unified 
treatment of the plasma process that is accurate at all temperatures and densities and also 
allows efficient numerical calculations. We introduce a simple analytic approximation to the 
dispersion equations for photons and plasmons which becomes exact in the classical limit, 
the degenerate limit, and the relativistic limit, and interpolates smoothly between these 
limits. Within the same approximation, we obtain simple expressions for the transverse, 
longitudinal, and axial vector components of the neutrino emissivity from the plasma process. 
The derivation of the analytic dispersion equations for photons and plasmons is presented 
in Appendix A. The effective neutrino-photon interaction that is responsible for the plasma 
process is discussed in Appendix B and the decay rate of a photon or plasmon into a neutrino 
pair is calculated in Appendix C. 

Because a plasma contains mobile charged particles, an electromagnetic wave propa- 
gating through the plasma consists of coherent vibrations of both the electromagnetic field 
and the density of charged particles. These coherent vibrations behave qualitatively differ- 
ently from electromagnetic waves in the vacuum in that there are longitudinal waves as well 
as transverse waves, and they propagate at less than the speed of light. The quantization of 
the electromagnetic waves in a plasma gives rise to a spin-1 particle with 1 longitudinal and 
2 transverse spin polarizations. It is common in the literature to refer to all 3 polarization 
states as "plasmons" , to emphasize that their dispersion relations depend on the properties 
of the plasma. The longitudinal and transverse modes are then awkwardly labelled "longi- 
tudinal plasmons" and "transverse plasmons". While the longitudinal mode owes its very 
existence to the plasma, the transverse mode simply has its dispersion relation at low fre- 
quencies modified by the plasma. For this reason and for the sake of concise terminology, the 
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label "plasmon" will be reserved in this paper for the longitudinal mode, and the transverse 
modes will be called "photons" whether they are propagating in a plasma or in the vacuum. 

The dispersion relations for photons and for plasmons depend on the temperature 
T and the net density n e of electrons minus positrons. The general expression for the net 
electron density as a function of T and the electron chemical potential p is 

n e (T,fi) = ^ J™ dpp 2 (n F {E) - np(E)^j . EB (1) 

Throughout this paper, we use units in which h — c — ks — 1- The integral in ([I]) is over the 
momentum p of electrons or positrons, E = \Jp 2 + ml is their energy, and m e is the electron 
mass. The Fermi distributions for electrons and positrons respectively are 

ME) = eiE .,l T+1 , (2) 

Given the net electron density n e , the chemical potential is determined by inverting ([]]) to 
obtain p(T, n e ) as a function of T and n e . 

In the vacuum, photons cannot decay into neutrino pairs because they are massless. 
Their dispersion relation is u 2 = k 2 , and the phase space available for the decay is propor- 
tional to uj 2 — k 2 = 0. The qualitative behavior of the dispersion relations in a plasma is 
shown in Figure 1. The upper and lower solid curves are the dispersion relations w t {k) for 
photons and ue(k) for plasmons at a temperature of T = 10 11 K = 8.6 MeV and an electron 
density corresponding to pj\i e = 10 12 g/cm 3 . (The quantity p/fi e is the mass density of 
protons in the plasma and is related to the net electron density n e by p/p e — m p n e , where 
m p is the proton mass.) As k — > 0, the dispersion relations u t (k) and uj(_{k) both approach 
the plasma frequency u p , which is given by 

"2 = v So dp p ~e{ 1 ~ \ y2 ) { nF(E) + nF{E) ) ' (4) 
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where v = p/E is the velocity of the electrons or positrons. At large k, the behavior of the 
photon dispersion relation is u>t{k) 2 — > A; 2 + m 2 , where m 4 is the transverse photon mass: 

ml = ^ £°dp£(n F (E) + n F (E)^j . (5) 

Comparison with (|j) reveals that the transverse mass lies in the range uj p < m t < \J<5/2 cu p . 

As /c increases, the dispersion relation Ut(k) for plasmons eventually crosses the light cone 
uj = k at a point fc max given by 

It satisfies u; p < k max < oo and represents the maximum momentum for which a plasmon can 
propagate. As discussed in Appendix A. 2, the expressions and @ for cu^, m 4 , and 

fc^a are correct to first order in the electromagnetic fine structure constant a ~ 1/137.036. 
Since the photon dispersion relation satisfies u t (k) > k for all k and the plasmon dispersion 
relation satisfies u>e(k) > k for k < k max , both photons and plasmons can decay into neutrino 
pairs. 

The production rate of neutrino pairs from the plasma process is sensitive to the 
precise form of the photon and plasmon dispersion relations. In 1961, Tsytovich wrote 
down general integral equations for the dispersion relations which include the effects of 
electrons in the plasma to first order in a. In their pioneering work on the plasma process 
in 1963, Adams, Ruderman, and Woo QT]] repeated the general integral equations, but they 
used the following simple dispersion relations in their numerical work: 

u t (k) 2 = u 2 p + k 2 , < k < oo , (7) 

uj e (k) 2 = uj 2 p , < k < u p . (8) 

At k — 0, these dispersion relations have the correct value u p at all temperatures and 
electron densities. At nonzero k, they are accurate only at temperatures and densities where 
the electrons are nonrelativistic. In 1967, Baudet, Petrosian, and Salpeter |4| improved on 



3 



the dispersion relations (|7|) and (pi) by including the first relativistic correction: 

u 2 k 2 

u 2 = u 2 p + k 2 + -± — , < k < oo , (9) 
5 u t 

2 3 cu? A; 2 
p 

where to>i is given by the integral 



uf = ^+^^ } 0<k<Ju 2 + 3u 2 /5, (10) 

O Up 



co 2 = ^ jT rf P £ - (n F{E ) + n F (E)^j . (11) 

It lies in the range < U\ < w p . At small k, the dispersion relations @ and (0) have the 
correct behavior to order k 2 at all temperatures and densities, but they are not accurate at 
large k if electrons are relativistic. 

In the subsequent 24 years, the dispersion relations (0) for photons and (|T0| ) for 
plasmons were used in all calculations of the plasma process || ^| It was pointed out 
by Braaten in 1991 that they could lead to significant errors at temperatures or densities 
where electrons become relativistic. In the relativistic limit, the correct dispersion relations 
are the solutions to the following transcendental equations: 

, , 9 9 3to> 2 / uf — k 2 u t , u t + k\ , , , 

9 9 3u; 2 f up , up + k \ , . , 

^ = ^TFU k «ZJ3i- 1 j' 0£,I<00 - (13) 

These dispersion relations were first derived by Silin in 1960 using kinetic theory ||, and 
were rederived in 1982 by Klimov and by Weldon using thermal field theory methods 0. 
Expanding the right sides of ( |i~2"|) and (|T3| ) in powers of k and using the fact that Ui = u p 
in the relativistic limit, one finds that they agree with (El) and ( [TP]) to order k 2 . They 
differ significantly at large k. For example, the transverse mass in the relativistic limit is 



m t = J 3/2 u p , while @ gives the value J 6/5 u p . Also, the maximum plasmon momentum 



approaches infinity in the relativistic limit, while ([K]) gives a value of J8/5 u p for k 



In a subsequent paper by Itoh et al. [jnj , the energy loss from the plasma process was 
calculated for a relativistic plasma using the dispersion relations for a degenerate plasma at 
zero temperature. The use of the zero-temperature dispersion relations is valid only if T is 
negligible compared to the electron Fermi energy. Since this condition is not satisfied by the 
hottest stellar plasmas, the calculation of Ref. |I(J can be improved by properly taking into 
account the effects of temperature on the dispersion relations. The general equations for the 
dispersion relations at zero temperature were obtained in closed form by Jancovici in 1962 
1 1 ij and used in the calculations of Ref. [10], but they are too lengthy to reproduce here. 
Jancovici also gave simplified equations for the dispersion relations: 



10+ 



k 2 + u 2 



p 2v 2 F k 2 



u? — Vpk 2 u t ujt + Vpk" 
log 



2v F k u>t — v F k 



< k < oo 



(14) 



u f 



3ujf ( uj£ ui + vpk 

log r - 1 



p v F k 2 \2v F k uj£ — vpk 



< k < k r , 



(15) 



where v F = pp/E F is the Fermi velocity, E F is the Fermi energy, and p F = \J E F — m 2 is 
the Fermi momentum. The maximum plasmon momentum is 



kr, 



3 / 1 . l + v F 
— U— log 1 

v F \2vp 1 — v F 



■nl/2 



LO r , 



(16) 



The dispersion equations ( P^ ) and ( p!5|) have been rederived using thermal field theory meth- 
ods [T2[]. In Ref. [IT|, these simple dispersion equations were only claimed to be valid for 
k « p F and uj « E F . However, as shown in Appendix A. 2, they are in fact valid for all 
k. They take into account correctly all effects of electrons in the plasma to first order in a. 
Thus they could have been used in the calculations of Ref. QTUJ without any loss of accuracy. 
Expanding the right side of flT4]) and flTSP for small k and using the fact that uj\ = v F uo p in 
the degenerate limit, we see that they agree with (P) and fllOD to order k 2 . In the relativistic 
limit where v F — > 1, ( |14|) and ( fTED reduce to the relativistic dispersion relations (0) and 
(PI). 



Aside from the numerical complications of solving integral equations, the general 
dispersion relations for photons and plasmons given in Ref. || also suffer from a technical 



difficulty in that uj t (k) and u>e(k) become complex- valued if the temperature or the Fermi 
energy is sufficiently high. The imaginary parts appear when uj p exceeds the threshold 2m e for 
decay of a photon into an electron-positron pair in the vacuum. This threshold is unphysical, 
because the rest mass of an electron in a relativistic plasma is significantly greater than in the 
vacuum. The effects of the plasma on the electron and photon dispersion relations are such 
that the decay of a photon or plasmon into an e + e~ pair is always forbidden by energy and 
momentum conservation |y|. The problem of complex- valued dispersion relations persists 
in the degenerate limit, although the threshold is increased to uo > Ep + m e , corresponding 
to the production of a positron at rest and an electron at the Fermi surface. The effects of 
the unphysical process 7 — > e + e~ may be numerically small, but it would be preferable to 
have dispersion relations from which they are absent altogether. 

In Appendix A. 2, we derive dispersion relations for photons and plasmons which 
receive no contributions from the unphysical process 7 — > e + e~. This is achieved without any 
sacrifice in accuracy: all effects from electrons and positrons in the plasma are still included 
to first order in a. The resulting expressions (§), and (||) for the plasma frequency, 
the transverse photon mass, and the maximum plasmon momentum are much simpler than, 
but just as accurate as, the corresponding expressions that follow directly from the general 
integral equations. The dispersion relations remain real-valued at all temperatures and 
electron densities. Their behavior at small k is consistent with (|9|) and (|I0D, they reduce to 
( |T2"D and ( |IBD in the relativistic limit, and they reduce to ([14]) and ( |T5| ) in the degenerate 
limit. 

The dispersion equations given in Appendix A. 2 are integral equations involving 1- 
dimensional integrals over the momenta of electrons and positrons. Only in the classical 
limit, the degenerate limit, and the relativistic limit can the integrals be evaluated ana- 
lytically. At intermediate temperatures and electron densities, calculating the dispersion 
relation u t (k) or u>e(k) requires finding, for each value of k, the zero in uj of an u;-dependent 
integral. The numerical solution of such a complicated dispersion equation is too inefficient 
for many applications, such as calculating the energy loss from the plasma process. We 
have therefore developed an approximation to these dispersion relations that is remarkably 
accurate at all temperatures and electron densities, but is simple enough to be used for prac- 
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tical calculations. The derivation of these dispersion relations is presented in Appendix A. 4. 
They are correct to order a in the classical limit, the degenerate limit, and the relativistic 
limit, and they provide a smooth interpolation to intermediate temperatures and electron 
densities. They are also correct to order k 2 at small k for all temperatures and electron 
densities. 

Our approximate dispersion relations are as simple as the degenerate dispersion rela- 
tions given in (|TJ) and (fr5|). Given the plasma frequency (^) and the freqency lo\ given by 
(P"TD, we define a parameter v*\ 



(17) 



It lies in the range < v* < 1, and has the intuitive interpretation of a typical velocity of 
electrons in the plasma. Our approximate dispersion relations uj t (k) and ue(k) are obtained 
by solving the following equations which depend on v*: 



0J< 



k 2 + u. 



p 2v 2 k 2 



co 2 -v 2 k 2 L0 t , U) t + vjz 

1 n ^^lOg 



2v*k 



< k < oo , (18) 



UJ, 



2 3u;f ( uj£ 1 ^ uj£ + v*k 



log : 



p v 2 k 2 \2v*k cue — v*k 
The maximum plasmon momentum is 



< k < k r , 



(19) 



3 / 1 . 1 + ^* . 

~5 I n l°g " 1 



2v* 



1 — v* 



1/2 



(20) 



The transverse photon mass is 



m t 



3 

2^! 



2v, 



-log 



1 + 



1/2 



P ■ 



(21) 



These expressions satisfy 2m 2 + (1 — v^k^^ = 3u 2 . The dispersion equations flUf ) and (jlQf) 
are correct to order a for all k in three limiting cases that are studied in Appendix A. 3. In 

the classical limit, flTT| ) gives t>* = ^5T/m e . The dispersion equations ( |T8"D and (|I9"D, when 

expanded to first order in agree with (|9|) and (p~0|) . In the relativistic limit, v* = 1 and 



( |IBD and fll9|) reduce to the relativistic dispersion equations ( |T2D and (|T3|). In the degenerate 
limit, is equal to the Fermi velocity up and the dispersion equations reduce to ( |I~4"D and 
(15). Our approximate dispersion relations are also correct to order k 2 at all temperatures 



and electron densities. Expanding the right side of (|T9|) and (|18|) in powers of k, we find that 
to order k 2 they reduce to @ and fllPl). At large k, the dispersion equations (|l8j) and (|I9| ) are 



correct to order a only in the classical, degenerate, and relativistic limits. For example, the 
expressions fl20|) and fl2l|) for A; ma:E and m t are not identically equal to the general expressions 



(^|) and (||), but the differences are found empirically to be remarkably small. 

To calculate the energy loss from the plasma process, we require the photon and 
plasmon dispersion relations uj t {k) and Ue(k), the corresponding residue factors Z t {k) and 
Zg(k), and the axial polarization function ILa(u, k) evaluated at the photon dispersion re- 
lation uj = uJt(k). General expressions for the residue factors and for n^(u;,fc) are given 
in Appendices A and B. The use of these general expressions in calculations such as the 
energy loss is very cumbersome. Fortunately, the methods used to derive the approximate 
dispersion relations ( |i~8"D and (plj) can also be used to derive compact analytic expressions 
for Z 4 , Z/, and H(w t , k) that are correct to order a in the classical limit, the degenerate limit, 
and the relativistic limit and are also correct for small k at all temperatures and electron 
densities. As shown in Apppendix A. 4, the transverse and longitudinal residue factors can 
be approximated by 

z (k) = 2 ^ 2 K 2 -^ 2 ) m 

A ) " 3o£w? + K + k 2 )(uj 2 - vlk 2 ) - 2u 2 (u 2 - k 2 ) ' [ZZ) 

Z (k) - 2 ^ ~ V * k ^ (23) 

As shown in Appendix B, the axial polarization function H(w, k) evaluated at u = Ut(k) can 
be approximated by 
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1-k 2 3^-2K 2 -P) 



U A (u; t , k) = u A k , |2 _ ai2 , 2 — ^ ~ ■ (24) 
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Its behavior as k — * is n^(o; t , A;) — > u A k, where the coefficient is 

u A = ^ jHdp i (l - (n F (£) - n F (E)^j . (25) 

The emissivity Q of a plasma is the rate of energy loss per unit volume. To calculate 
the emissivity from the plasma process, one must first calculate the rates T t (k) and T^(fc) 
for the decay of a photon and a plasmon of momentum k into a vv pair. These rates are 
calculated in Appendix C under the assumption that the number density of neutrinos and 
antineutrinos remains negligible. The emissivity is then obtained by integrating over the 
phase space of the photon or plasmon, weighted by the number density and by the energy. 
Summing over the 3 polarization states of photons and plasmons and over the neutrino types 
v = u e , u T , the total emissivity is 

Q = E/ ( 2 nsMk)) ujt(k) T t (k) + n B (u e (k)) coe(k) . (26) 

The number densities of photons and plasmons are given by the Bose distribution: 

= • (27) 

The total emissivity Q can be separated into the vector (Qt) and axial vector (Q A ) compo- 
nents of the photon contribution and the plasmon contribution (Ql)'- 

E C v j I dk k 2 Z t (k) (u t {kf - k 2 ) n B (u t (k)) , (28) 

^^jg^fe/o dkk2z t( k )("t(k) 2 -k 2 ) U A (cu t (k),k) 2 n B (cu t (k)) , (29) 
Ql = (j2 C v) jf"" dk k " M k ) ^{k) 2 (^{k? - k 2 f n B (u e (k)) . (30) 



The coefficients Cy and C A depend on the neutrino type and are given in Appendix B. The 
combinations that arise in the plasma process are 

ECy = t _ 2sin 2 ^ + 12 sin 4 6 W « 0.911, (31) 
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4 



(32) 



The formulas for Qt and Qa were first correctly given up to the factor of C v in Ref. [14 
The correct values of Cy were first included in Ref. ||. The formula (|30"D for Qt differs 
from that given in Refs. H and [K| because the definition of the longitudinal residue factor 
Zt(k) differs by a factor of {ujf — k 2 )/uo 2 . The formula ( p9[ ) for Qa supercedes that given in 
Ref. pT5|| , in which the factor Z t (k) was omitted and the photon was assumed to satisfy the 
simple dispersion relation (|7]). 



The important momentum scales in the emissivity integrals appearing in (£g), fl29|), 
and fl30D are the plasma frequency u p (which enters into the dispersion relations, the residue 
factors, and rU(cut, k)), and the temperature T (which occurs in the Bose distribution). The 
integrals can be simplified in the limiting cases T » uj p and uj p « T. We first consider 
the high temperature limit T » uj p . The integral (p8|) for the transverse emissivity is 
dominated by momenta k of order T. Since T » uj p , the photon dispersion relation can be 
approximated by uo 2 = k 2 + m 2 . We can therefore replace the factors of uj t (k) 2 — k 2 in ( |28| ) 
by m 2 , and set ut = k everywhere else in the integrand. The integral over k can then be 
evaluated analytically in terms of the Riemann zeta function: (3 ~ 1.202057. The result is 



4C 3 m 6 t T 3 



(33) 



The integral ( p9|) for the axial vector emissivity is also dominated by momenta k of order T. 
Applying similar approximations as were used for the transverse emissivity, it can also with 
some effort be evaluated analytically: 



Qa - (Eci)^ 



1 , l + v, 
log 1 



2v 



1 



i 2 



2T 

m t 



(34) 



Up to logarithmic factors, the axial vector emissivity is suppressed relative to Qt by a factor 
of uj\/T 2 . In the limit T » u p , the integral (p0|) for the longitudinal emissivity is dominated 
by momenta k of order uj p . This is obvious in the nonrelativistic limit, because the upper 
limit k max is of order uj p . In the relativistic limit where k max approaches infinity, the factors 
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Zi(k) and (uji(k) 2 — k 2 ) 2 both approach zero very rapidly for k » u p , providing a cutoff on 
the integral of order u p . If k is of order lu p , then ue(k) is also of order ui p . Using u>e(k) <<T, 
the Bose factor in (^) can be simplified to ne(u^) — > T/uj£. The integral can still not be 
calculated analytically, but by dimensional analysis it must be proportional to uj^ T. The 
longitudinal emissivity then has the form 

(^)^A(v,)^T. (35) 

The coefficient varies slowly from 8/105 in the nonrelativistic limit (v* — > 0) to 0.349 

in the relativistic limit (u* — > 1) 0. The longitudinal emissivity is suppressed relative to the 
transverse emissivity by a factor of u 2 /T 2 . 



We next consider the low temperature limit T « uo p . The integrals in ( JZ8"D and ( |29"| ) 
are dominated by momenta k « uj p . The Bose distribution can therefore be approximated 
by a Gaussian in k : n B {u t {k)) -> e- <JJ f/ T ea;p(-^ / (0)fc 2 /2T), where <(0) = (1 + ^/5)/^. 
Everywhere else in the integrands, we can set u t = u p and ignore k relative to u p . The 
integrals can then be evaluated analytically, with the results 

«-(? c ^*( i+ ¥) ^ < 12 Tm ^ rlT • < 36 ' 

The axial vector emissivity is suppressed relative to Qt by a factor of uP^TjuJt. A similar 
approximation can be applied to the longitudinal emissivity provided that T << v 2 uj p . 
The Bose distribution in ([50) can then be approximated by a Gaussian in k: ns{uJe{k)) — > 
e -^p/ T exp{-u'l{0)k 2 /2T), where ^'(0) = (3/5)v*/u p . The expression for the emissivity then 
reduces to 

«i - (E Cl) ^ f- (§*) ^ ^ ^ . (38, 

At relativistic electron densities (t>* — > 1), the longitudinal emissivity is smaller than the 
transverse emissivity only by a factor of y2. 
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In Figure 2, the components Qt, Qa, and Ql of the emissivity are shown as a function 
of the proton mass density p/fi e — m p n e at a temperature T = 10 11 K = 8.6 MeV. Also 
shown as dotted lines are the corresponding emissivities calculate using the O-temperature 
dispersion relations as in Ref. flQfl . While they give the same results at high densities, 
there are significant discrepancies at the lowest densities considered in Figure 2. The shapes 
of the dispersion relations at T = and T = lO 11 -^ are very similar, since the electrons 
are relativistic over the entire range of densities shown. The discrepancies in Figure 2 
therefore arise primarily from the difference in the value of the plasma frequency. At the 
highest densities shown, uo p is determined primarily by the Fermi energy, so the O-temperature 
dispersion relations provide a good approximation. At the lowest densities shown, u p is 
determined primarily by the temperature. The emissivities Qt and Ql therefore become 
almost independent of the density, as can be seen in Figure 2. 

The formulas for the dispersion relations and the neutrino emissivity that have been 
presented above were derived for a plasma of electrons and positrons in a uniform positively 
charged background that cancels the net charge of the electrons and positrons. In a stellar 
plasma, the cancelling charge is actually provided by protons and heavier ions. These charged 
particles will also give contributions both to the electromagnetic polarization functions and 
to the effective photon-neutrino interaction that add to those of the electrons and positrons. 
The effects of protons and heavier ions are negligible for stellar plasmas as long as they are 
nonrelativistic. At the highest densities shown in Figure 2, the protons are relativistic and 

their effects are significant. At T = 10 11 K and p/fj, e = 10 14 g/cm 3 , electrons and protons 
are degenerate with both having a Fermi momentum of 238 MeV. Their contributions to 
ujp are proportional to their Fermi velocities, which is 1 for electrons and 0.246 for protons. 

Thus the protons increase the plasma frequency by about 12%. The resulting effect on the 
emissivity will be much larger, since it scales like a large power of the plasma frequency. 
The modifications to the formulas for the emissivities that are required in order to take into 
account the effects of protons are given in Appendix C. 

We have developed a unified approach for calculating the neutrino energy loss from 
the plasma process at all temperatures and densities. We have introduced compact equations 
for the dispersion relations for photons and plasmons that are correct to order a for all k in 
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the classical, degenerate, and relativistic limits and are also correct to order k 2 at small k 
for all temperatures and electron densities. Compact expressions were also obtained for the 
other quantities that are required to calculate the transverse, longitudinal, and axial vector 
components of the neutrino emissivity. Most previous numerical studies of the energy loss 
from the plasma process have concentrated for simplicity on the case where the neutrino 
density remains negligible. In supernova explosions, neutrinos become trapped inside a 
neutrinosphere and the neutrino emissivity can be suppressed by Pauli blocking effects. Our 
approach should allow the efficient numerical investigation of such effects. It would also 
be useful to extend our approach to other neutrino emission processes, to axion emission, 
and to other particle physics proceeses that can play an important role in astrophysics over 
enormous ranges of temperatures and densities. The investigation of all of these processes 
would benefit from a unified treatment that remains accurate even under extreme variations 
of temperature and density. 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, under Grant DE-FG02-91-ER40684. 



13 



A Photon and plasmon dispersion relations 



A.l Dispersion relations to order a 

The effects of a plasma on the propagation of photons and plasmons is determined by the 
electromagnetic polarization tensor W(K). If interactions with electrons and positrons 
are taken into account to leading order in the electromagnetic coupling constant a, the 
polarization tensor is 



U^(K) = 16™ J ^ ^f(E) + n F {E) 



P ■ K(P»K U + K*P V ) - K 2 P»P V - (P ■ Kfg^ 
X {P-K) 2 -{K 2 ) 2 /A ' (39) 

where = (u>, k), P^ = (E,p), K 2 = uj 2 — k 2 , and P ■ K = Eu —p - k. The integral is over 
the momentum p of electrons and positrons with energy E = ^p 2 + m 2 . The polarization 
tensor satisfies T^IP" (K) = 0, which is a consequence of gauge invariance. The transverse 
and longitudinal polarization functions Pitiyj, k) and Pie(u, k) are 

Pi t (u, k) = ^ (S ij - k l k j ) U ij (uj, k) , (40) 

Pi t {u,k) = U 00 {uj,k) . (41) 

They are related to the standard transverse and longitudinal dielectric functions e t and en 
by e t = 1 — Pi t /uj 2 and q = 1 — Pi e /k 2 . 

In order to construct the effective propagator D^ v {oj, k) for the electromagnetic field, 
it is necessary to choose a gauge. The most convenient choice for treating plasma effects is 

— * — * 

the Coulomb gauge defined by V ■ A = 0. In Coulomb gauge, the nonzero components of 
the effective propagator are 



^■^ F-Fi,!^) ' (42) 
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D ij {u,k) 



8 l3 -k l k J ) . 



(43) 



oo 2 — k 2 — Pit(uJ, k) 

The dispersion relations w t {k) for photons and u>i(k) for plasmons are the locations of the 
poles in the effective propagator: 



n oo/ n ^(k) 2 ^(*0 

D™(co,k) -> — r — — as cj->u^(/c), 



D ij (cu,k) 



k 2 u 2 — LUe(k) 2 
Z t (k) 



lo 2 - uj t {kf 



(5 ij - k r k j ) as u -> w t (fc) 



(44) 
(45) 



The locations of the poles are independent of the choice of gauge [ITS]. 



The equations ( PHI ) and (FJSp also define the residue functions Z^(fc) and Z t (k), which 
determine the strength with which a plasmon or photon of momentum k couples to an 
electromagnetic current: 



Z t {k) 



1 - ^±{uj t {k),k) 



du 2 



Z e (k) 



k 2 



uj t (kf 



dPi e 
duo 2 



-i -l 



(w t (k),k) 



(46) 



(47) 



The longitudinal residue Z^{k) defined by (|47| ) differs by a factor of (u 2 — k 2 )/u 2 from 
that which was used in Refs. || and ||10|| . The residue of a pole in u 2 of D^ v {uj,k) can 

be identified as e M (fc)e^(/c)*, where e^(k) is the polarization 4-vector for the appropriate 
propagating mode. These modes are conveniently labelled by the helicity A, which is the 

component of the angular momentum in the direction of A;: A = for plasmons, A = ±1 for 
photons. From (|44] ) and (f4~5|) , we identify the polarization 4- vectors to be 

u t (k) 



A = 0) = -^V^)^ )" 
e"(£,A = ±l) = ^Uk) (0,e±(k 



(48) 
(49) 



where e+(k) and e~(k) are orthogonal to k and normalized so that e±(k) ■ e±{k)* = 1. The 
polarization 4-vectors (Pi8"D and (|9|) satisfy the Coulomb gauge constraint 



k ■ e{k, A) = . 
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(50) 



They are used in Appendix C to calculate the decay rate of a photon or plasmon into neutrino 
pairs. 

The dispersion relations 0J t (k) for photons and 0Ji{k) for plasmons are the solutions 
to the equations 

u t {kf = k 2 + Pi t (cu t (k),k) , (51) 

Mkf = ^^Pi e (u; e (k),k). (52) 

If the complete order-a expression ( |39D for the polarization tensor is used to calculate these 
dispersion relations, they become complex-valued when the temperature or electron density 
is large enough that the plasma frequency u p exceeds 2m e . The imaginary part is propor- 
tional to the decay rate for a photon or plasmon into an electron-positron pair: 7 — > e + e~ . 
These imaginary parts are unphysical, because the plasma effects that give the photon a 
nontrivial dispersion relation also change the energy-momentum relation for the electron. 
These corrections are such that the decay 7 — > e + e~ is always forbidden by energy and 
momentum conservation [|13|| . The unphysical effects from this forbidden process also reveal 
themselves in the analytic behavior of the real parts of the dispersion relations. Thus they 
can not be eliminated simply by defining the dispersion relations oj t {k) and ue(k) to be the 
real parts of the poles in the effective propagator. 

A. 2 Removing effects of electron pair production 

To eliminate the unphysical effects of the forbidden decay 7 — > e + e~, we exploit the fact 
that there is a separation of momentum scales between the particles whose propagation is 
significantly modified by the plasma and those that are responsible for the plasma corrections. 
For example, in the relativistic limit, the dominant contributions to the electromagnetic 
polarization tensor in fl39| ) come from electrons and positrons with momentum p of order the 
temperature T or the Fermi momentum pp, whichever is larger. However plasma corrections 
to the electron propagator are significant only for momenta that are smaller by a factor of 
\f(x. In the nonrelativistic limit, plasma corrections to the electron propagator are always 

negligible. This justifies the use of the vacuum energy- momentum relation E = Jp 2 + m 2 
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in calculating IP 1 '. 

The separation of momentum scales also allows us to simplify the expression (|39| ) 
for the polarization tensor by dropping the term (K 2 ) 2 /4 in the denominator. The physical 
interpretation is that this corresponds to calculating plasma corrections using the forward 
scattering amplitudes for electrons and positrons in the vacuum [|I7|]. The mathematical 
approximation that is required is 

\u 2 - k 2 \ « 2E\u-v-k\ , (53) 

where v = p/E is the velocity of the electron or positron. We first discuss the case where k is 
of order uj p . Since the dispersion relations u t and ue are then also of order u p , the inequality 
( |53"|) requires uj p « E. This is always satisfied in the nonrelativistic limit where E is of 
order m e . In the relativistic limit, the inequality uj p « E fails to be satisfied only in a region 
of the electron or positron phase space whose contribution to the integral in (|39D is down 
by a factor of w 2 /T 2 in the high temperature limit and a factor of Lo^jE 2 F in the degenerate 
limit. In either case, the contribution is suppressed by a factor of a. We next discuss the 
case in which k is much greater than uj p . The solutions to the dispersion relations always 
have uj 2 — k 2 of order u 2 . Except in the relativistic limit, the inequality ([53]) then reduces 
to uj 2 /k << E, which is even more easily satisfied than the condition uo p « E discussed 
above. In the relativistic limit, the inequality ( |53"D may fail for electrons with velocity v 
within an angle of order u p /k of the momentum k, but this region of phase space gives a 
contribution to the integral in fl39|) that is suppressed by at least a factor of a. Thus the 
effects of the (K 2 ) 2 term in the denominator of ( |3"5| ) are always at most comparable to the 
order a 2 corrections to the polarization tensor. It can therefore be dropped without any loss 
in accuracy. 

Dropping the (K 2 ) 2 term in the denominator of (|3^) eliminates the effects of the 
unphysical process 7 — > e + e~, so that the dispersion relations for photons and plasmons 
remain real-valued at all temperatures and densities. It also results in far simpler expressions 



for the transverse and longitudinal polarization functions defined in ([|0) and (41). The 
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angular integrals can be evaluated analytically, and the polarization functions reduce to 



4a r™ p 2 (u 2 u 2 - A; 2 w 

*) = ^ r dp ^ f > g ^ - 1 - 4=*) Ue) + M ^ . (55) 

7r jo i? \ffc uj — vk uj 1 — v 2 k 2 J \ J 

Inserting ([54]) and ( |55|) into (|5T|) and (J52|) and solving these equations numerically, one finds 
that the dispersion relation uj t {k) remains real- valued for all k and u>e(k) remains real- valued 
for those values of k that satisfy u>i(k) > k, corresponding to the timelike propagation of 
a plasmon. As k — > 0, both dispersion relations approach the plasma frequency, which is 
given by @). At large k, the behavior of the photon dispersion relation is governed by the 
transverse photon mass given by (||). The plasmon dispersion relation crosses the line uj = k 
at the point k max given by @. 



nj,(£) + n F {E) , (54) 



A. 3 Limiting cases 

The polarization functions ([54]) and (|55D can be evaluated analytically in three limits: the 
classical limit, the degenerate limit, and the relativistic limit. This allows the equations 
( |5"ID and ( p^ ) for the dispersion relations to be written in closed form. We also obtain 
analytic expressions for the plasma frequency (f|) and the residue functions ( ^6|) and (f47|). 
For completeness, we also give in each case the solution of ([I]) for the chemical potential /i 
function of the net electron density n e . 

Classical limit. In the classical limit, the plasma is nonrelativistic (T << m e ) and nonde- 
generate (m e — \x » T). The Fermi distribution for electrons can be approximated by the 
Boltzmann distribution e^~ E ^ T , and contributions from positrons can be ignored. The net 
electron density, including the first correction proportional to T/m e , is 

n e (T,v) = e^V T ^=(m e Tf 2 U + ~) ■ (56) 

V In 6 \ o m e / 

This can be solved trivially for the chemical potential n(T, n e ) as a function of the temper- 
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ature and the net electron density. The plasma frequency oo p is 



Using (pq) to eliminate \x in favor of the net electron density, the expression (|57|) reduces to 



47ra;n e / 5 T 

p m e V 2m t 

The polarization functions, including the first correction proportional to T/m e , are 

/ k 2 T \ 

U t (u,k) = ul 1 + — — , (59) 

/ p k 4 T \ 

Ut( Ul k) = u? v l- 2 + 3^— . (60) 



The resulting dispersion equations are 



k 1 T 



k l + ui 1 + — , < fc < oo , (61) 
'V wfmj 



/ k 2 T \ i 

uf = 1 + 3— — , < fc < ^1 + 3T/m e w p . (62) 



These are identical to the dispersion equations @ and (0) with uj\ = JhT/m e lu p . The 



transverse mass for the photon is m t = y 1 + T/m e uo p . The residue factors are 

^(fc) = i - -At— , (63) 

W) = 4 f 1 - 6 4-) • (64) 
u 2 \ u 2 m e J 

Degenerate limit. The degenerate limit is the limit of low temperature T « 
jj, — m e . In the limit T = 0, The Fermi distribution for electrons reduces to a step function: 
rip(E) = 1 for E < fi and np(E) = for E > \i. The net electron density is 

n e (T = 0,/i) = ^ (fi 2 -m 2 e ) 3/2 . (65) 
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The chemical potential at T = is called the Fermi energy: fi = Ep. The solution to (|65| ) 



for the Fermi momentum pp = J Ep — as a function of n e is 



Pf 



(37r 2 n e ) 



1/3 



(66) 



The plasma frequency u p reduces in the T = limit to 



4a 
3^ 



U P = —Pf v f , 



(67) 



where the Fermi velocity is Vp = Pf/ Ep. The polarization functions are 



H t (uj,k) = u), 



3u 2 (_ u 2 — vpk 2 uj , u + vpk" 



p 2v 2 F k 2 



uj 2 2vpk ^u — vpk) 



(6£ 



Ut(u,k) = u 2 — 

^ «p \2vpk uj — vpk 



u . u + Vpk ^ N 
log r - 1 



The dispersion equations reduce to ( |i~4"D and (|i~5|) . The transverse photon mass is 



(69) 



2^ 



1-vj 
2v F 



log 



1 + uf 
1 - v F/ 



1/2 



(70) 



The inverses of the residue factors are 



-i 



3 ui 



2v 2 F k 2 \2 



3u 2 — Vpk 2 u t ^ u t + vpk y 
2u\ 2vpk LUt — vpk 



(71) 



z - x = 3 u p f uj 

2 vpk 2 \uj — Vpk 2 2vpk^ b ut — vpk) 



U£ ui + vpk 
log- 



(72) 



The dispersion equations (|i~4l) and (|i~5f ) can be used to eliminate the logarithms from (|71~1) 
and ([F2|). The resulting algebraic expressions for Z t and Z\ have the form (|22]) and (p3l), 
with = t/p. 

Relativistic limit. The relativistic limit is the limit of either high temperature T » m e 
or high density fi >> m e . In the limit m e = 0, the net electron density is 



n e (T,/j,) = 7^^(V + ^T 2 ) ■ 



(73) 
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The solution for the chemical potential as a function of n e is 



(jl(T, n e ) 



\ 



1/3 



\ 



1/3 



where = (37r 2 n e ) 1//3 . The plasma frequency is 



(74) 



2 4(3; / 2 1 2rri2 

= 3^ + 3* T 



(75) 



The polarization functions are 



U t (u,k) = u; 2 



3u^ 
2k 2 



^ ur — Ar u; ^ + A; 
a; 2 2A; g a; — k 



(76) 



3w - lot 



a; + k 



2k uj — k 



- 1 



(77) 



The dispersion equations reduce to ([12]) and (|i~3|). The maximum plasmon momentum is 



infinite and the transverse mass is m t = v/3/2 dj p . The inverses of the residue factors are 



^r 1 



1 3 uJp /3 3cu 2 - A; 2 + k 

~ 2 1 2 V 2 2w 2 2fc ° g uj t - k 



(78) 



Z7 1 



3wi 



^7 



0^ , ^ + fc 

2P lu;?-ifc 2 2A? ° g ^ - A; 



(79) 



The dispersion equations (|T2|) and (|13D can be used to eliminate the logarithms from (]78f) 
and (|79|) . The resulting algebraic expressions for Z t and have the form (|22|) and ( p3|) with 



u* = 1. 



A. 4 Analytic Approximation 

The dispersion relations (|5l| ) and ( 52]) are integral equations and solving them numerically 
can be computationally time-intensive. It is therefore desirable to have analytic approxima- 
tions to these equations. To derive such equations, we begin from the expressions ( ]54|) and 
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fl55|) for the polarization functions and integrate by parts: 



„ . lX 4a f 00 , p 3 { 1 feu 2 uj 2 - v 2 k 2 uj , cu + f/A] d ( , ' 



(80) 

P^,k) = --1 \- 2 fag -1 - „,(S) + H,(B) ■ (81) 



If the integrals over the momentum p of the electrons and positrons were dominated by a 
single velocity t>*, then the factors in the curly brackets could be evaluated at v = v* and 
pulled outside the integral. This is in fact the case in the classical limit, the relativistic limit, 
and the degenerate limit. In the classical limit, all the electrons have velocities v near 0. In 
the relativistic limit, electrons and positrons all have velocities v* = 1. In the degenerate 
limit, the factor dnpjdp is sharply peaked at the Fermi velocity vp. After pulling the 
expressions in curly brackets out of the integrands in ( |5C| ) and Q8"T|), the remaining integrals 
over p are proportional to the square of plasma frequency. This can seen by using integration 
by parts on M): 



~ r d p^( n ^ + • (82) 



dp \ 

The resulting analytic expressions for the polarization functions are 



2 3 ( uj 2 oj 2 — v 2 k 2 to , to + v*k s 



P,^,k) = ^ ( log—j-ll. (84) 



Inserting these expressions into fl5T|) and (|52|) , we obtain the dispersion equations flPED and 



For suitable choices of the velocity u*, the expressions (|83 ) and (|8~4" ) will be accurate 
in the classical, degenerate, and relativistic limits. The parameter v* can also be chosen so 
that they are correct at small k for all temperatures and electron densities. At small k, the 
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general expressions fl54|) and (||IJ) for the polarization functions reduce to 



j-j)^\(n F (E) + n F (E)) , (85) 



«,(„,*) « ^/; dP ^{(i_^) £ + £}(„, (B) + MB) ) . (86, 

Expanding out the analytic expression ( |S"3]) for the transverse polarization function in powers 
of k, we find that it agrees with (|85"D to order k 2 provided that we take v* = u)\/u p , where u)\ is 
given in (|TTD . We find also that the analytic expression ( pi]) for the longitudinal polarization 
function agrees with (|5!5| ) to order k 4 for the same value of v*. 

We can also obtain compact analytic expressions for the residue factors Z t and Zf. 
Inserting (|8~3"D and (|4]) into the formulas (^) and (|47|) , we obtain 

; 2 /3 gt^-igg 
2 y2fc 2 \ 2 2cu 2 2^£; ~° w t - 



7 -i -i 3^/3 3a; 2 - vjk 2 u t u t + v*k 
Z t = 1 " ol^fe o JO oTTT S ~ — ) , (87) 



3 ^ 2 / a; 2 ^ ^ + ^fc \ 

' 2^ 2 V^!-^ 2 2^ S u t -v.k) " 1 j 

The dispersion equations (|18|) and (|19D can be used to eliminate the logarithms from (|8 



and ( [38]) , resulting in the simple algebraic expressions (p2) and (|23|). 



B Effective photon-neutrino interaction 

The decay of a photon or plasmon into neutrino pairs proceeds through an effective photon- 
neutrino interaction. This effective interaction arises from the electromagnetic coupling of a 
photon to electrons or positrons in the plasma, together with the weak interaction coupling 
of the electron or positron to a neutrino pair. It can be summarized by an effective vertex 
r aM for the interaction of the photon field with the neutrino current uj a (l — 75)^: 

r-M) = ^^=(c vW) (1,^(1,^ 



+ 9 C 



C v Pi t (u, k) (5 ij - titi) + C A TI a (lu, k) (ie l]m k m )} g j »j , (89) 
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where = (oJ,k) is the 4-momentum of the photon. This effective vertex satisfies the 
identity T^K^ = 0, which guarantees the gauge invariance of the interaction. The functions 
Pit{uJ, k) and Pii(uj, k) in (|89"D are the transverse and longitudinal electromagnetic polariza- 
tion functions defined in fl40|) and (41). The axial polarization function calculated to leading 



order in a is 

, , N K 2 r d 3 p 1 / , s . \ P-Kuj-K 2 E , . 

U A (.,k) = j ^ -{n F[ E) - ME)j {p , K)2 _ (K2m , (90) 

where K 2 = uj 2 — k 2 and P ■ K = Eu — p- k. The coefficients Cy and Ca depend implicitly 
on the neutrino type. The vector coefficients are 

C v = 2sm 2 9 w + ~ for v e , (91) 

Cy = 2sin 2 9 w - - for v^, v T , (92) 
where 9\y is the weak mixing angle: sin 2 9\y ~ 0.226. The axial vector coefficients are 

C A = \ for v e , (93) 

C A = -- for v T . (94) 

The axial polarization function (|90| ) agrees with the expression used in Ref. [[01 , provided 
that the nonrelativistic dispersion relation (^) is used to set K 2 = uj 2 . 

For frequencies uj greater than 2m e , the expression (^) for the axial polarization 
function has an imaginary part that arises from the production of e + e~ pairs. As discussed 
in Appendix A, this unphysical behavior can be eliminated without any loss of accuracy 
by dropping the term (K 2 ) 2 /4 in the denominator of (|90|). The resulting expression for 
HA{w,k), when evaluated on the photon dispersion relation u = u> t (k), is still correct to 
leading order in a. Having dropped the (K 2 ) 2 /A term in the denominator of (|90|) , the axial 
polarization function reduces to 

„ . , N 2a u 2 - k 2 f°° , p 2 ( uj , u + vk uj 2 - k 2 \ ( 

(95) 
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Its behavior for small k is Ua{u, k) — > u>Ak, where uoa is given in (|25|). 

The expression ( f?5|) for the axial polarization function can be evaluated analytically 
in the classical limit, the degenerate limit, and the relativistic limit: 

Classical limit. At leading order in T/m e , the axial polarization function reduces to 

W) = e<— y^Jl a (H) 1 " M^fl . (96) 

V 7T \ ITLf I UJ 



Eliminating the chemical potential and using the expression ([57]) for the plasma frequency 
in the classical limit, fl9"6|) simplifies further to 

„ , 1N wl k(oo 2 -k 2 ) , . 

Degenerate limit. At T = 0, the axial polarization function reduces to 

U A (uJ,k) = —p F - — -log --1 . (98) 

7r k \2vpk uj — Vpk J 



After setting tu = u> t (k), the logarithm can be eliminated using the dispersion relation (14). 
The resulting expression for U.j^(u t , k) is (plj) with v* = vp and uja = (2a/37r)pFV F . 

Relativistic limit. For m e = 0, the axial polarization function reduces to 

rr < 7^ 2a ^ 2 - k 2 ( u u) + k \ 



After setting lu = u t (k), the logarithm can be eliminated using the dispersion relation (|12] 



The resulting expression for II^u^, k) is (|24j) with v * = 1 and c^a = (2a/37r)yU. 



We can derive an analytic approximation to the axial polarization function using the 
same methods that were applied to the transverse and longitudinal polarization functions in 
Appendix A. Starting with the expression (|95|) and integrating by parts, we obtain 



„ . ,. 2a J 2 - k 2 r°° , p 3 [ 1 / uj , + vk , , • , , 



(100) 
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In the classical, degenerate, and relativistic limits, the integral is dominated by a single 
momentum v*. The factor in curly brackets can therefore be evaluated at v* and pulled 
outside the integral. The remaining integral over p is proportional to uja, which after applying 
integration by parts to the expression in (j25l) can be written in the form 



lot r°° t? A ( \ 

"* = -5f/. * h * H £) - ME) ■ (101) 



Our final expression for the axial polarization function is 

rr i n uj 2 - k 2 3 f u uj + v*k \ 

This expression is not only correct for all k in the classical, degenerate, and relativistic 
limits, but it has the correct behavior for small k at all temperatures and electron densities. 
The axial polarization function H^(u t , k) evaluated at the photon dispersion relation can be 
further simplified by using the dispersion relation ([19|) to eliminate the logarithm in (|102|) , 



resulting in the analytic expression given in (|24|). 



C Decay rate of photon or plasmon 

In this Appendix, we calculate the rates T t (k) and Ti(k) for a photon or plasmon of momen- 
tum k to decay into neutrino pairs. We begin by writing down the matrix element M. for 
the decay into a neutrino and antineutrino of momenta p\ and p2- 

M = ^ (T^A A )) fi (Pi)7«(l - 75MP2) , (103) 

where e^(k, A) is the polarization 4-vector given in (|49|) for photons or dig) for plasmons, 
T afl is the effective photon-neutrino vertex given in (p9|), and u(pi) and v (p^) are the spinors 
for the neutrino and antineutrino. The matrix element must be squared and integrated over 
the phase space of the neutrino and antineutrino: 

- ^!w?kl^k^ mp > +F *- K)m2 - (104) 
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The subscript A on T and on u is either t, corresponding to helicity A = ±1, or £, correspond- 
ing to helicity A = 0. In the delta function, Pi = (pi,pi), Pi = (^2,^2), and K = (ux^k), k) 
are the 4-momenta of the neutrino, antineutrino, and photon or plasmon, respectively. In the 
expression ( |104j ), we have assumed for simplicity that the number density of neutrinos and 
antineutrinos is negligible. Otherwise the phase space integrals in (|104|) must be weighted 
by appropriate Pauli blocking factors. The only dependence in the matrix element (|103|) on 
the momenta pi and p 2 is in the spinor factor. After multiplying the spinor factor by its 
complex conjugate, it can be expressed in the form of a Lorentz tensor: 

M(j0i)7 a (l - 75)w(P2)v(p 2 )7 /3 (l - 75)m(pi) 
= 8 {if Pi + P 2 °Pf - Pi • P 2 g a P - ie a ^P ltl P 2u ) . (105) 

In the absence of Pauli blocking factors, the integral over the phase space of the neutrino 
and antineutrino in ( |104[ ) can be carried out analytically, with the result 



' f ( '' l>2 1 (2n) 4 5\P 1 + P 2 - K) «(pi) 7 Q (l- J 5 )v(p 2 ) - ->b)«(ft) 



(2vr) 3 2piy (2vr) 3 2p 2 

= ±-{K a K p - K 2 g al3 ) , (106) 

where K 2 = w\ — k 2 . Since the effective vertex T a ^ satisfies the identity K a T afl = 0, only 

the K 2 g af3 term in (|106|) contributes to the decay rate ( |104| ). The decay rate then reduces 
to 

Uk) = - ^ Ux{ * ] * { - k2 (r*%(M)) (r ap e»(k,\)Y . (107) 

To complete the calculation of the decay rate, the effective vertex T afJ- in (^) must 

be contracted with the appropriate polarization 4- vector e^(k, A) and the energy uj must be 
evaluated at the corresponding dispersion relation u) t {k) or u)i{k). For the plasmon, only 
the first term in (|89|) contributes. By the plasmon dispersion equation fl52|), Hi(w, k) can be 



replaced by k and the contraction of T a ^ with the polarization vector (Pf ) reduces to 

r°*{u>e(k),k) e„(£,0) = C v "^== \/~Ze(k) k uj e (k) • (10* 
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For the photon, it is the last term in the effective vertex fl39| ) that contributes. By the photon 
dispersion equation (|5l|) , Ht(uj, k) can be replaced by uj t (k) 2 — k 2 and the contraction of T afJ - 
with the polarization vector (^) reduces to 

Gr 



r^(ut(k),k) eii (k,±i) = -^J^icv^kf-k^^Uki 

- C A U A (uj t (k), k) (o, ik x e ± (k)Y ) . (109) 



Squaring the expression ( pLU9| ) and ( |108| ) and inserting into ( |107| ), the decay rates of the 
photon and plasmon into neutrino pairs reduce to 

r t (k) = z t (k) ^ ~ e [c 2 v (u t (kf - k 2 ) 2 + c\ u A {u t {k\k)^ . (no) 
r <(*) = °v tUt z ^ k ) M^) 2 ~ k2 T ■ ( m ) 



The expressions ( |111|) and (|110|) for the decay rates of photons and plasmons include 



only the effects of electrons and positrons in the plasma. It is straightforward to include also 
the effects of protons provided that the plasma frequency remains small compared to 700 
MeV, so that form factor effects can be neglected. The effects of protons on the dispersion 
relations and on the effective neutrino-photon vertex must both be included. The expressions 
( |110|) and (|111|) for the decay rates are replaced by 



r*(*) = ^ r z t ( ^—^-( (c v p#\u t ,k) - h v P$\uH,kj) : 



+ (C A r4 e) (u; t , k) - h A IT? K k)) 2 ) ' (112) 



Uk) = -^f- Z, Ue ^ k2? (C v P#\ Ui , k) - h v PiPfa, k)) 2 , (113) 

where Pi\ , P%f \ and IT^ are the contributions to the transverse, longitudinal, and axial 
polarization functions from electrons, while Ptf\ Pif , and are the corresponding 
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contributions from protons. The proton terms have the same form as the electron terms 
(|83D , fl84|), and ( |102| ), except that the parameters lu p , v* = uii/ujp, and to a are calculated 
using the integrals @, (|TT1), and (|2~5|) with the electron mass m e replaced by the proton 
mass m p . The coefficients hy and Iia that describe the interactions of the proton with 
neutrinos are 

- 2 9w + - 
2 



hy = —2 sin 2 Q w + - for v P , v„ , z/,. 



e j "/i > t j 



(114) 



hy 



1^9 A for v e , v^, v T , 



;ii5) 



where qa ~ 1.26. The dispersion equations (|5lD and (|52]) must also be modified to include 
the effects of protons: 



= k 2 + Pti e) (u t ,k) + P£?\u u k) 



;ii6) 



k 2 



(Pif\u e ,k) + Pif{u>t,k) 



(117) 



The residue factors Z t in ( |112| ) and in ( |113|) are then given by the expressions fl37|) and 
fl38|), except that in addition to the electron term on the right side, there is also a proton 
term of the same form but with appropriate values for the parameters u p and v*. 
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Figure Captions 



1. Dispersion relations u(k) for photons (upper solid curve) and plasmons (lower solid 
curve) at temperature T = 10 11 K and proton mass density p//x e = 10 12 g/cm 3 . 

2. Transverse (T), longitudinal (L), and axial vector (A) components of the neutrino 
emissivity (in units of erg/s/cm 3 ) as a function of the proton mass density p/fi e ( m 
units of g/cm 3 ) at the temperature T = 10 11 K (solid curves). Also shown are the cor- 
responding emissivities calculated with the O-temperature dispersion relations (dashed 
curves) . 
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